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We discuss the implications for gravitational wave detectors of a class of modified gravity theories 
which dispense with the need for dark matter. These models, which are known as Dark Matter 
Emulators, have the property that weak gravitational waves couple to the metric that would follow 
from general relativity without dark matter whereas ordinary particles couple to a combination of 
the metric and other fields which reproduces the result of general relativity with dark matter. We 
show that there is an appreciable difference in the Shapiro delays of gravitational waves and photons 
or neutrinos from the same source, with the gravitational waves always arriving first. We compute 
the expected time lags for GRB 070201, for SN 1987a, and for Sco-Xl. We estimate the probable 
error by taking account of the uncertainty in position, and by using three different dark matter 
profiles. 
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I. INTRODUCTION 

The direct detection of gravitational waves from astro- 
physical sources would enable us to open a new window 
into the universe and get insights which are complemen- 
tary to electromagnetic astronomy Many ground- 
based interferometric detectors such as LIGO, VIRGO, 
GEO600 and TAMA have been online for several years. 
In October 2007, LIGO completed a long science run to 
collect one year of coincident data at design sensitivity [2] 
and the VIRGO detector also joined this science run in 
the last five months. During the latest LIGO science run, 
the sensitivity of the 4 km Hanford and Livingston LIGO 
detectors to detect binary neutron-star coalescence with 
mass 1.4 Mq with signal to noise ratio greater than 8 
(after averaging over all orientations and sky positions) 
was about 15 Mpc Analysis of the latest LIGO and 
VIRGO data for gravitational wave (GW) searches from 
a wide variety of sources is in progress Q . 

An important science goal pursued is the search for im- 
pulsive transient GW signals from sources with electro- 
magnetic and/or neutrino counterparts. Some examples 
of such sources include core-collapse supernovae, gamma- 
ray bursts (GRBs), soft gamma-ray repeaters (SGRs), 
pulsar glitches, low mass X-ray binaries, blazar flares, 
optical transients, etc Q. These "triggered" searches al- 
low us to get better sensitivity for a given false alarm 
rate as compared to an all-sky search at all times and 
to design custom-made analysis algorithms taking into 
account our knowledge of the source astrophysics. Con- 
versely, there has been a proposal to look for optical and 
infrared counterparts at the time of coincident GW burst 
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candidates [5|. An overview and benefits of such triggered 
searches carried out by the current interferometric grav- 
itational wave detectors are reviewed in Ref. There 
have been proposals to determine neutrino mass using 
simultaneous neutrino and GW observations from core- 
collapse supernova [6(. Similar triggered GW searches 
will also be important for the future LISA experiment 0] . 

In all present and past triggered searches for gravi- 
tational waves, the analysis is done by looking at the 
data from GW detectors within a narrow time window 
(of about hundreds of seconds) around the time of the 
electromagnetic trigger. With this assumption, one can 
detect gravitational waves only if the propagation time 
of photons/neutrinos is the same as that of gravitational 
waves. In general relativity, photons, neutrinos and grav- 
itational waves propagate on the same null geodesies. 
Hence the total time of propagation is the light travel 
time delay plus the Shapiro time-delay due to intervening 
matter 8]. For electromagnetic waves, Shapiro delay has 
been detected in a wide variety of systems such as radar 
ranging to Venus, Doppler tracking of Cassini spacecraft 
and in binary pulsars [9| . From the relative arrival times 
of photons and neutrinos from SN 1987a, we also know 
that the Shapiro time-delay for neutrinos is the same as 
that for photons to within 0.5% p^-Hlj. 

The conventional view is that general relativity de- 
scribes gravity on cosmic scales. If this is so, the gravi- 
tation of stars and gas is not sufRcient to account for the 
velocity dispersions in clusters 11211 , or for the rotation 
curves of spiral galaxies Il3l.ll4l. ll5l. or for the weak lens- 
ing in galactic clusters [l6lll7l.[lg^ll9l. [20l . [2]| . Big Bang 
Nucleosynthesis severely limits the extent to which the 
deficit can be made up of unseen but ordinary matter [2^ ; 
the remainder must consist of an exotic, nonrelativistic 
substance which has never been detected except gravita- 
tionally. This dark matter must vastly predominate over 
ordinary matter. For example, only about one fifth of our 
galaxy's mass is made up of normal matter, with the rest 
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being composed of dark matter [23] . Thus the dominant 
contribution to the Shapiro delay for photons from GRBs 
and other sources is due to the gravitational potential of 
the intervening dark matter. 

None of the proposed dark matter candidates have 
been detected, either directly in a laboratory exper- 
iment or indirectly through their annihilation prod- 
ucts [13, HE [13 ■ This prompts the suspicion that 
perhaps it is gravity which must be modified, rather 
than the universe's inventory of nonrelativistic matter. 
Of course that would invalidate the assumption which is 
the basis for all current and proposed GW searches from 
sources seen in photons and neutrinos [1, [^ . A previous 
study [H, [1^ has considered the consequences for grav- 
itational wave detection of a certain class of modified 
gravity theories known as Dark Matter Emulators. In 
this paper we correct a mistake in the original work [1^ 
that led to the wrong sign for the effect, and we work out 
explicit results for three interesting sources. 

Section II defines and motivates Dark Matter Emula- 
tors. In Section III we review three popular dark mat- 
ter profiles which these models are designed to obviate. 
Section IV computes the expected time lag between the 
arrival of the pulse of gravitational waves from some cos- 
mic event and its optical or neutrino counterpart. In 
Section V we give explicit results for three sources of in- 
terest. Section VI gives a very brief discussion of other 
alternate gravity models, and our conclusions comprise 
Section VII. 



II. DARK MATTER EMULATORS 

Certain regularities in cosmic structures suggest mod- 
ified gravity. One of these is the TuUy-Fisher relation, 
which states that the luminosity of a spiral galaxy is 
proportional to the fourth power of the peak velocity in 
its rotation curve [30| . If luminous matter is insignifi- 
cant compared to dark matter, why should such a rela- 
tion exist? Another regularity is Milgrom's Law, which 
states that the need for dark matter occurs at gravita- 
tional accelerations of oq — 10""'^'' m/s^ [3l|. A third 
regularity is that uq also seems to give the internal accel- 
erations of pressure-supported objects ranging over six 
orders of magnitude in size — from massive molecular 
clouds within our own galaxy to X-ray clusters of galax- 
ies [13. 

A modification of Newtonian gravity which explains 
these regularities was proposed by Milgrom in 1983 ^3^. 
His model. Modified Newtonian Dynamics (MOND), was 
soon given a Lagrangian formulation in which conserva- 
tion of energy, 3-momcntum and angular momentum are 
manifest 34] ■ However, there was for years no success- 
ful relativistic generalization which could be employed 
to study cosmological evolution. Even in the context of 
static, spherically symmetric geometries. 



the early formulation of MOND fixed only B{r), not 
A{r). It was therefore incapable of making definitive pre- 
dictions about gravitational lensing. 

A relativistic extension of MOND has recently been 
proposed by Bekenstein [35]. This model is known as 
TeVeS for "Tensor- Vector-Scalar." In addition to repro- 
ducing the MOND force law at low accelerations, TeVeS 
has acceptable post Newtonian parameters, and it gives 
a plausible amount of gravitational lensing [s^. When 
TeVeS is used in place of general relativity -I- dark mat- 
ter to study cosmological evolution, the results are in 
better agreement with data than many thought possi- 
^40|. The model does have problems 

times 



ble [361 |37|,|3&,p...,p.^j. .... ^^^^ ..^.^ 

nth stability [4l|, [44] . The Bullet Cluster is sometir 



cited as a fatal blow for the model [45] but opinion on 
this differs ^4^, HE [1^ , and this system in any case poses 
problems for dark matter [H, \^ . 

What concerns us here is the curious property of TeVeS 
that small amplitude gravitational waves are governed, 
as in general relativity, by the metric ji^j/, whereas mat- 
ter couples to a "disformally transformed" metric which 
involves the vector and scalar fields. 



(2) 



-B{r)c^dt^ + A{r)dr^ + r^dVt^ 



(1) 



The Scalar- Vector- Tensor gravity (SVTG) theory pro- 
posed by Moffat also has different metrics for matter and 
small amplitude gravitational waves [H, [4^. The ap- 
pearance of this feature in two very different models is 
the result of trying to reconcile solar system tests with 
modified gravity at ultra-low accelerations. Solar system 
tests strongly predispose the Lagrangian to possess an 
Einstein-Hilbert term [H^, [5lj . On the other hand, failed 
attempts to generalize MOND ^53] have led to a theorem 
that one cannot get sufficient weak lensing from a stable, 
covariant and purely metric theory which reproduces the 
TuUy-Fisher relation without dark matter |53j . Hence the 
MOND force must be carried by some other field, and it 
is a combination of this other field and the metric which 
determines the geodesies for ordinary matter. However, 
the dynamics of small amplitude gravitational waves are 
still set by the linearized Einstein equation. This simple 
observation makes for a sensitive and generic test. 

We define a Dark Matter Emulator as any modified 
gravity theory for which: 

1. Ordinary matter couples to the metric 5^1^ that 
would be produced by general relativity -|- dark 
matter; and 

2. Small amplitude gravitational waves couple to the 
metric g^^ produced by general relativity without 
dark matter. 

Now consider a cosmic event such as a supernova which 
emits simultaneous pulses of gravitational waves and ei- 
ther neutrinos or photons. If physics is described by 
a dark matter emulator then the pulse of gravitational 
waves will reach us on a lightlike geodesic of g^^, whereas 
neutrinos and photons travel along a lightlike geodesic of 
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If significant propagation occurs over regions that 
would be dark matter dominated in general relativity 
then there will be a measurable lag between arrival times. 

Currently the only observational constraint on the 
speed Vg of gravity relative to that of ordinary mat- 
ter Win derives from the consequences of gravitational 
Cherenkov radiation from particles moving faster than 
gravity [H, [s^. From observations of the highest en- 
ergy cosmic rays Moore and Nelson infer the bound, 
Vm — Vg < 2 X 10~^^c Although the original study of 
Dark Matter Emulators [1^ in fact violated this bound, 
that was the result of incorrectly choosing the dimen- 
sional constant in a certain logarithm. In the next section 
we show that the speed of gravity is always greater than 
that of light for Dark Matter Emulators. A discussion 
of the Shapiro delay calculation in some other alternate 
gravity theories can be found in Refs. [sgI. [57|. 



Profile 




ro 


Tc 


Isothermal 


3.98 days 


4.00 kpc 


219 kpc 


NFW 


60.8 days 


16.7 kpc 


N.A. 


Moore 


51.8 days 


29.5 kpc 


N.A. 



TABLE I: Dark Matter Profile Parameters for the Milky Way 
Galaxy from Ascasibar et al. [60| ]. 



difference in propagation times between the two geome- 
tries! We will henceforth consider p{r) to be the density 
of dark matter. 

The potentials AA{r) and AB{r) can be given simple 
expressions in terms of the mass function. 



III. THREE DARK MATTER PROFILES 

We shall specialize to static, spherically symmetric dis- 
tributions of dark matter, consistent with the invariant 
element ([T|). It is well to bear in mind that hierarchical 
structure formation will not necessarily result in spher- 
ically symmetric distributions [58|. There is even evi- 
dence that the dark matter halo conjectured to surround 
our own galaxy is not spherical (soj . 

For a pressureless, static, spherically symmetric system 
the Einstein equations take the form. 
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rA V / 



Vb 



A-l 



SttG 
& 
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(3) 



If the potential B(r) goes to a constant at infinity we can 
choose the time units so that equation ^ has the exact 
solution, 



B[r) ~ exp 



dr 



(5) 



However, our study requires only small corrections to 
A{r) and B(r), 



A(t) = 1-1- A^(r) 



B{r) = 1 + AB(r) . (6) 



This not only simplifies ((3][4|), it also means we can dis- 
pense with the contribution of ordinary matter to the 
mass density p(r). The reason is that we are computing 
the difference in propagation times along null geodesies 
between the same points in two different geometries. The 
geometry felt by gravitational waves is sourced only by 
ordinary matter, while the geometry felt by photons and 
neutrinos is sourced by the sum of ordinary matter and 
dark matter. At first order in the mass density, the ef- 
fect of ordinary matter cancels out when computing the 



M(r) = 47r / dr' p{r') 
Jo 



The linearized solution of ([3]) is. 



(7) 



AAir) 



8nG 



\ , ,2 , 2G M{r) 
dr r^p{r ) = — 



Note that AA{r) is positive semi-definite. From 
find AB{r), 



(8) 



we 



(4) AB{r) 



dr 



, AA{r') 



-AA[r) 



2G 

72 



dr' 



M{r') 



(9) 

Note that AB{r) is negative semi-definite and in fact less 
than or equal to —AA{r). This guarantees that gravita- 
tional waves travel faster than photons or neutrinos, so 
there is no problem with the bound of Moore and Nel- 
son [s^ . 

Our study requires the dark matter density functions 
for our own galaxy and (for the most distant source) 
for the Andromeda galaxy. We took these in the form 
of fits to three popular density profiles whose analytic 
forms are presented at the end of this section. Given the 
current rough quality of the observational data, a Dark 
Matter Emulator that reproduced the potentials AA{r) 
and AB{r) for any of these profiles would be judged suc- 
cessful. One can therefore regard the slightly different 
time delays that result as one measure of the theoretical 
uncertainty. The fits for our own galaxy appear in Ta- 
bleland were taken from the study by Ascasibar, Jean, 
Boehm and Knodlseder of the positron annihilation line 
from the galactic center The fits for Andromeda 

appear in Table |TT] and were done by Tempel, Tam and 
Tenjes 
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Profile 


87rGporg/c3 


ro 


Tc 


Isothermal 


1.88 days 


1.47 kpc 


117 kpc 


NFW 


48.6 days 


12.5 kpc 


N.A. 


Moore 


45.8 days 


25.0 kpc 


N.A. 



values of po and rg given by Ascasibar, Jean, Boehm and 
Knodlseder [g^l, along with = 219 kpc. This choice 
for Tc causes the ratio of the total masses of the Milky 
Way and Andromeda galaxies for the isothermal profile 
to agree with that of the NFW profile considered in the 
next subsection. 



TABLE II: Dark Matter Profile Parameters for the An- 
dromeda Galaxy from Tempel et al. [gH - 



B. NFW Profile 



A. Isothermal Halo Profile 

We shall use a variant of the Isothermal Halo Profile in 
which the density vanishes after a cutoff radius Vc [62...63j. 
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(10) 



Such a cutoff is inevitable, even in MOND, owing to the 
presence of other galaxies. Of course it is also necessary 
to make the potential Ai3(r) vanish at infinity. 
For r < Tc the mass function and potentials are, 



M{t) = A^p^rl I tan-i ( - ) - -— ^ 



Ayl(r) 
AB(r) 
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'l-^tan-i 
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ro> 3(r2+r2) 



(11) 
(12) 
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(13) 



For r > Tc the mass is constant and the (equal and op- 
posite) potentials fall off like 1/r, 



M(r)^\'KpQrl\ 3^ -tan 
A^(r) 



"(-) 



Jo ^tq, 
^Gparl \r,{2rl + Srg) 
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■ tan 
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(15) 
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'0^ 



-^tan-Yi:^) .(16) 
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It should be noted that Ascasibar, Jean, Boehm and 
Knodlseder used a simpler version of the isothermal pro- 
file without a cutoff radius re- This would cause the 
potential AB{r) to eventually become positive, which 
violates the bound of Moore and Nelson It also 

doesn't make any sense when one considers the effect 
of other galaxies. Because the isothermal profile is the 
most closely related to MOND we considered it impor- 
tant to include results for this profile, so we used the 



The NFW profile was the result of studying the equi- 
librium density profiles of dark matter halos in numerical 
simulations of structure formation l64i. 
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The associated mass function and potentials are, 

M(r) 

A^(r) 
AB{r) 
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(18) 

(19) 
(20) 



C. Moore Profile 

A later effort along the same lines showed a better fit 
to a density function which is more sharply peaked at the 
center [63, 
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The associated mass function and potentials are. 
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IV. THE SHAPIRO DELAY FOR GIVEN M{r) 

It is more convenient to convert the spatial coordinates 
from spherical (r, 9, (f) to Cartesian x', 



sm W cos 0, sm f sm 0, cos ( 



(25) 



The invariant element ([T|) has a simple expression in 
terms of these coordinates, 



c^dr +dx-dx- AB c'^dr + AA{f- dxY . (26) 



Before moving on we digress to note that specializing 
to ds^ — and identifying the velocity of photons and 
neutrinos as = dx/dt results in an equation for the 
speed of effectively massless, ordinary matter, 







-(1 + AB)C^ + ■ Vm + AAif ■ VmY 



(27) 



Now recall that we are ignoring the role of ordinary mat- 
ter in both the metrics of gravity (where it is the only 
part of the mass density) and of ordinary matter. It fol- 
lows that the speed of gravity is Vg is c. Treating to first 
order in the potentials and assuming AA{r) > and 
AB{r) < 0, we see that Wi„ — Wg < 0. 

We can express the invariant element ([26|) as the flat 
space contribution plus a perturbation, 



ds^ = {rj^i, + h^i,)dx^'dx'' 



(28) 



Comparing and allows us to read off the 3-1-1 
decomposition of the graviton field h^^^, 

hao = -AB , ho^^O and h^j ^ AAf'f^ . (29) 

One advantage of Cartesian coordinates is that the affine 
connection vanishes for the flat background. It is easy to 
give the first correction. 



pa 



per 



(30) 
(31) 



We need the null geodesic x^(t) that connects the 
spacetime points x^ = (0,a?i) and X2 = (ci,a?2). It obeys 
the geodesic equation. 



r + rv(x(T))x''r = 0, 

subject to the conditions. 



1 ! 



X^(0) = X' 



(32) 



(33) 
(34) 
(35) 



Of course we solve this perturbatively in the potentials. 
The zeroth order solution is of course the flat space result, 



X^{t)^x'I + Ax>^t . 



(36) 



Here the temporal and spatial components of the interval 
Ax'' are. 



Ax*' = \\X2 — Xi 



and Ax' = X2 — x\ 



(37) 



The first order corrections to the spatial components 
of the geodesic are, 

x\{t) J^dr'il- T')Ar%^(x + AxT^AxPAx" 

- dr' (r - r') Ar^p^ (xi + Axt) Ax" Ax'" . (38) 

Of course it is from the first order temporal correction 
that we infer the time lag. This correction is more com- 
plicated, 



x;(r) 



2Ax 



V (2^1) Ax'' Ax'" 



T 

Ax 



j dr' (1 - r')Ax' Ar%^ (x + Axr) Ax" Ax"" 
- dr' (t - r') Ar% (xi + Axr) Ax" Ax'" . (39) 

Ignoring ordinary matter makes the graviton geodesies 
identical to Xo {''')■ Hence the time lag between the ar- 
rival of gravitational waves and the arrival of photons or 
neutrinos is (to first order in the potentials), 

cAt EE xi{l) = I dr h^^ (xi + Axt) Ax'' Ax" . 

2Ax Jo ^ ' 

(40) 

Expression can be simplified a great deal further. 
First expand out the potentials. 



hp^Axf^Ax" = -AB Ax^ + AA{f- Aa 



(41) 



(Note that the time lag is positive semi-definite because 
AB < and AA > 0.) Now use relation jS]) for AB{r) 
in terms of AA{r) and partially integrate to reach the 
form. 



Ax . , 
cAt = -—AB{r2) 

-t-— / dTAA{r) 



T cir /r ■ Ax\ 2 
r dr \ Ax / 



Ax • xi . , , Ax ■ Xo . , , 



2 Ax ^ ' 2 Ax 



+Ax / drAAir) 1 



(xi • Ax)^ - rjAa 
Ax^r^ 



It is useful to define the constant C, 
1 



C: 



Ax2 



'5^Ax2 - (xi • Ax) 2 



Finally, we change variables from r to r, 

'Xi ■ Ax\ 2 



C2 . 



(42) 



(43) 



(44) 



(45) 
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Object 


R. Ascen. 


Decl. 


Dist. 


Milky Way 


17h45m 40s 


-29° 00' 28" 


7.94 kpc 


Andromeda 


OOh 42m 44s 


+41° 16' 09" 


778 kpc 


GRB 070201 


OOh 44m 32s 


+42° 14' 21" 


780 kpc 


SN 1987a 


05h 35m 28s 


-69° 16' 12" 


51.4 kpc 


Sco-Xl 


16hl9m55s 


-15° 38' 24" 


2.80 kpc 



TABLE III: Angular coordinates and distances for the Milky 
Way and Andromeda galaxies and for the three sources of this 
study. 



Assuming r2 < ri the result is, 

Ax ■ xi . , Ax ■ X2 . ^ 

cAt = — ABin) AS ra) 

2Ax ^ ' 2Ax ^ ' 

r\ 2GM{r) r /CAx\2 ^^^^ 

For ri < r2 we take the other root of the solution for 
T(r), which reverses the upper and lower limits in P5)l . 



V. RESULTS 

We have worked out explicit results for three typical 
sources at vastly different distances: GRB 070201, SN 
1987a and Sco-Xl. Their celestial coordinates are given 
in Table ITTll Table [ml also gives the centers of the Milky 
Way and Andromeda dark matter halos. 

GRB 070201 was a short hard gamma ray burst whose 
angular error box corresponded to a 0.124° quadrilateral 
which overlapped with the Andromeda galaxy (66j . Short 
hard gamma-ray bursts are believed to be caused by the 
mergers of two neutron stars or a neutron star and a black 
hole ^]. If GRB 070201 derived from such a merger, 
with masses close to 1.4 Mq and a reasonable orientation, 
the GW signal should have been seen if its distance was 
780 kpc [6^ . It is however possible that the GRB did not 
originate in the Andromeda galaxy, and in that case the 
signal from a compact object merger may be inaccessible 
to LIGO. 

No gravitational waves were found from a search done 
within ± 180 s time-window with the LIGO Hanford de- 
tectors around the time of this GRB 69. One in- 
terpretation of this null result is that GRB 070201 was 
a SGR flare [i^, [Tlj. However, it is also possible that 
physics is described by a Dark Matter Emulator, in which 
case the pulse of gravitational waves would have arrived 
long before the electromagnetic signal. Table IIVI gives 
our results for the time lag one would expect at the cen- 
tral position using each of the three dark matter density 
profiles. Although the time lags differ by as much as 69 
days, none of the lags is less than two years. 



Profile GRB 070201 SN 1987a Sco-Xl 

Isothermal 742 days 78.2 days 4.98 days 

NFW 804 days 74.8 days 4.88 days 

Moore 811 days 74.5 days 4.97 days 



TABLE IV: Time delays from three dark matter profiles for 
each of the three sources of this study. 

Table |V] considers another measure of the likely error 
by specializing to the isothermal profile and varying the 
angular position (at the fixed distance of 780 kpc) over 
the four vertices of the angular error box. In this case dis- 
tinct results are reported for the contributions from the 
Milky Way and Andromeda halos, which are of course 
independent at linearized order. As expected, varying 
the position has no effect on the contribution from the 
Milky Way halo but it can change the contribution from 
the Andromeda halo by as much as 15 days. We should 
however stress that if this delay calculation was done by 
assuming that this GRB is at a distance of 780 kpc. If 
this GRB did not originate in Andromeda, then the cal- 
culated delay would be much larger. 

SN 1987a was a core collapse supernova in the Large 
Magellanic Cloud at a distance of 51.4 kpc J^L Neu- 
trinos were observed by the Kamiokande-II [TJ, UM ^^'^ 
Irvine-Michigan-Brookhaven [tI, [tJ detectors. The op- 
tical signal arrived several hours later because photons 
must traverse the optically dense stellar environment [t^ ] . 
The total Shapiro delay for SN 1987A from the visible 
and dark matter distribution has also been independently 
estimated [13, [z3 to be between 0.29 to 0.36 years. If the 
oblateness of SN 1987a was in relation to that of the Sun, 
the current gravitational wave detectors would probably 
not have seen anything had they been operating at the 
time [tS^. However, advanced LIGO would detect such 
a supernova out to 0.8 Mpc [78]. This includes the An- 
dromeda galaxy, which doubles the expected rate and 
also ensures that the signal passes through dark matter 
dominated regions. Of course the effective coverage from 
neutrino detectors will remain limited to our galaxy and 
its satellites (zl, . 

Table IIVI gives our results for the expected time lag 
from a Dark Matter Emulator which reproduces each of 
the three dark matter profiles. These results include only 
the effect of the Milky Way halo. In contrast to the 
much more distant GRB 070201, the scatter between the 
various models for SN 1987a is much smaller — a mere 
2.7 days. 

Sco-Xl (located at a distance of 2.8 kpc) is one of the 
brightest Low Mass X-ray Binaries (LMXBs). LMXBs 
are potential sources of gravitational waves from r-modes 
gettin g e xcited due to accretion, or from a deformed 
crust |8ll . HI]. One proposed search is to look for 
coincidences between the data from LIGO and Rossi X- 
Ray Timing satellite Q. This search also assumes that 
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R. Ascension 


Dechnation 


AtMW 


AfMSl 


OOh 44m 32s 


42° 14' 21" 


407 dy 


335 dy 


OOh 46m 18s 


41° 56' 42" 


407 dy 


337 dy 


OOh 41m 51s 


42° 52' 08" 


407 dy 


322 dy 


OOh 42m 47s 


42° 31' 41" 


407 dy 


330 dy 


OOh 47m 14s 


41° 35' 35" 


407 dy 


338 dy 



TABLE V: Shapiro Delays for GRB 070201 from the Isother- 
mal Profiles of the Milky Way (AtMw) and Andromeda 
(AtM3i) at the central value of the angular position and at 
the four vertices of the error box. In all cases the distance to 
the burst was taken to be 780 kpc. 



gravitational waves and X-ray photons arrive at the same 
time. 

Table HVl reports the expected time lag for each of the 
three dark matter profiles, again from the Milky Way 
halo. Although the time lag is still easily observable at 
~ 4.9 days, the agreement between the three models for 
this source is excellent. The largest discrepancy is just 
.1 day. 

VI. OTHER MODIFIED GRAVITY THEORIES 

Since the 1970s, there have been various proposed tests 
of general relativity through gravitational wave observa- 
tions (See Ref. [SJI for a recent review). Most of these 
tests are in the strong field regime. In this section, we list 
some other non-GR gravity theories which also predict a 
non-zero time-delay between photons and GWs and are 
not yet ruled out through other observations. 

These are massive graviton theories and brane-world 
models. In massive graviton theories [ssj, the gravita- 
tional waves would arrive after the photons, with the 
delay being dependent on the graviton mass. However, 
the Moore-Nelson lower bound on the speed of gravita- 
tional waves imposes stringent constraints on the validity 
of massive graviton models. It should be noted that there 
is no consistent interacting theory for massive spin two 
particles which is limited to a finite number of fields [s^ . 

In various brane-world models, gravitational waves 
propagate faster than photons or neutrinos, depending 
on the curvature of the bulk [13, [H, 111 • Therefore, it is 
not possible to calculate model-independent time-delays 
for the three sources we considered in this paper. 

VII. CONCLUSIONS 

The power and generality of our analysis derives from 
ignoring the details of how a Dark Matter Emulator dis- 
penses with the need for galactic dark matter. We merely 
assume that it does, which implies that ordinary matter 



must couple to the metric that general relativity would 
predict with dark matter. The special characteristic of 
Dark Matter Emulators is that weak gravitational waves 
couple to the metric that general relativity would pre- 
dict without dark matter. Both of these metrics can be 
inferred from observation, and all geometrical quantities 
worked out, without regard to the details of specific mod- 
els. Although a Dark Matter Emulator is not the only 
conceivable way of evading the no-go theorem [5^ while 
preserving solar system tests js^l , it is the only way that 
has so far been given a concrete realization. 

If dark matter does not exist and the observed cosmic 
motions and lensing instead derive from a Dark Matter 
Emulator then the assumption upon which all triggered 
gravitational wave searches are based breaks down. In 
this case the optical or neutrino identification of a plau- 
sible gravitational wave source would not imply a simul- 
taneous pulse of gravitational waves but rather that such 
a pulse occurred earlier. Even for nearby sources such as 
Sco-Xl (at about 2.8 kpc) the gravitational waves would 
arrive almost five days earlier. For a source in the An- 
dromeda galaxy the time difference would be over two 
years. 

It is obviously premature to proclaim that the failure of 
triggered searches to reveal any coincident gravitational 
wave pulse implies that physics is described by a Dark 
Matter Emulator. But if plausible sources continue to 
produce null results this possibility has to be considered. 
In that case the key question becomes the accuracy with 
which one can estimate the expected time lag. Some 
measure of this is given by the spread in Table HVl for dif- 
ferent reasonable dark matter profiles. Table fVl considers 
variations in the angular position, and there will be com- 
parable results for varying the much less well-determined 
distances. Based on these analyses it seems unlikely that 
the uncertainty can be reduced below the level of a few 
percent. This has important implications for the way 
data needs to be kept and for the types of searches that 
should be contemplated. 

Of course a gravitational wave signal might be loud 
enough to show up in all-sky untriggered searches [90j . 
In that case the gravitational wave signal would trigger 
electromagnetic and neutrino searches. A single gold- 
plated event in which the counterpart signal arrived after 
a plausible delay would be powerful evidence in support 
of Dark Matter Emulators. Conversely, a single detec- 
tion of coincident signals would rule out the entire class 
of Dark Matter Emulators. This is a novel way of us- 
ing gravitational wave detectors to test alternate gravity 
theories in the ultra-weak field regime. Indeed, this is an 
ideal test because Dark Matter Emulators do not change 
aspects of the tensor component of a gravitational wave 
signal such as the number of polarizations or, to any rea- 
sonable accuracy, the travel time between Earth-bound 
detectors. So there need be no change in the data anal- 
ysis algorithms that would be used in any case. 
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